LEMNISCATES DO NOT SURVIVE LAPLACIAN 

GROWTH 



D. KHAVINSON, M. MINEEV-WEINSTEIN, M. PUTIN AR, 
AND R. TEODORESCU 



1. Introduction 

Many moving boundary processes in the plane, e.g., solidification, 
electrodeposition, viscous fingering, bacterial growth, etc., can be math- 
ematically modeled by the so-called Laplacian growth [9,13]. In a nut- 
shell, it can be described by the equation 

(1-1) V(z) = d n g m (z,(), 

where V is the normal component of the velocity of the boundary dQ(t) 

d 

of the moving domain C 1R 2 ~ C, z G d£l(t), t is time, — denotes 

on 

the normal derivative on d£l(t) and gn(t)(z,() is the Green function 
for the Laplace operator in the domain Q(t) with a unit source at the 
point ( G Q(t). Equation (11.11) can be elegantly rewritten as the area- 
preserving diffeomorphism 

(1.2) S(z t z e ) = l, 

where 3 denotes the imaginary part of a complex number, dQ(t) : = 
{z := z (t, 8)} is the moving boundary parametrized by w = e %e on the 
unit circle and the conformal mapping from, say, the exterior of the 
unit disk D + := {\w\ > 1} onto Q(t) with the normalization z(oo) = 
(,z'{oo) > 0. 

The equation (11.21) . named Laplacian growth or the Polubarinova - 
Galin equation in modern literature, was first derived by Polubarinova- 
Kochina [11] and Galin [7] in 1945, as a description of secondary oil 
recovery processes. 

This equation is known to be integrable [10], and as such possesses 
an infinte number of conserved quantities. More precisely, it admits 
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conserved moments c n = J z n dx dy, where n runs over either all non- 

Q(t) 

negative or all non-positive integers depending on whether domains 
Q(t) are finite or infinite. At the same time (11.21) admits an impressive 
number of closed-form solutions. 

For the background, history, generalizations, references, connections 
to the theory of quadrature domains and other branches of mathemat- 
ical physics we refer the reader to [4, 8-10, 12, 13] and the references 
therein. 

In section §2] of this paper, we show that any continuous chain of 
polynomial lemniscates of order n: T t : = {|P(z, £)| = 1}, P(z,t) = 

n 

a(t) Y\ ( z — <\j-(£)), where a(t) is real-valued, is destroyed instantly 

3=1 

under the Laplacian growth process described in fll.ip . with Q(t) = 
{\P(z,t)\ > 1}, ( = oo, unless n — 1, Ai(t) = const and {T t = dQ(t)} 
is simply a family of concentric circles. Here the roots Xj(t) of P(z,t) 
are all assumed to be inside Q' t := {\P(z,t)\ < 1}, so Q and Q' are 
simply connected. 

This result shows that unlike quadrature domains (cf. [4,12]) that are 
preserved under the Laplacian growth process, lemniscates for which all 
the roots of the defining polynomial are in f2£ are instantly destroyed, 
except for the trivial case of concentric circles. This, incidentally, agrees 
with a well-known fact — cf. [5] - that lemniscates which are also 
quadrature domains must be circles. The proof of the theorem for the 
case of Laplacian growth is given in §21 

In §3] we extend the result of §|2] to all the growth processes that are 
invariant under time-reversal and for which the boundary velocity is 
given by 

(1-3) V(z) = x(z)d n g m (z,(), 

with x{ z ) is a bounded, real, positive function on T t . Invariance under 
time-reversal is defined here in the following way: if the boundary T t+( it 
is the image of T t under a map f^ dt) : z t G Y t i-> z t+dt G T t+dt , then 

f{t+dt-dt) ° f(t,dt) = I- 

We conclude with a few remarks in §H 

2. Destruction of Lemniscates 

Theorem 2.1. Suppose that a family of moving boundaries T t , (where 
t > is time), produced by a Laplacian growth process, is a family of 

n 

polynomial lemniscates {\P(z,t)\ = 1}, where P(z,t) = a(t) Y[ [ z ~ ^j{t)]> 

3=1 
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and all Xj(t) are assumed to be inside T t . Then, n = 1 and Ai = const, 
i.e., T t is a family of concentric circles. 

Proof. Let Sl t = {z : \P(z,t)\ > 1}, D+ = {\w\ > 1}. The function 
(f(t) : Q t — > D + , w = (p(z,t) = ^/P(z, t), where we choose the branch 
for the n— th root so that (p'(t, oo) > 0, maps Qt conformally onto D + , 
(p(t, oo) = oo. It is useful to note that on T t , P(z, t) = w n , \w\ = 1 and 
does not depend on t. This is because for any two moments of time t, r, 
we have w(t)(.) = w(t) o n(t, r)(.), where k is a Mobius automorphism 
of the disk. In our case, n(t, t)(oo) = oo, so K(t,r)(z) = e ia z,a G R, 
but since it also fixes the argument at oo, k is the identity. 

Therefore, we have (where, as is customary, we denote the partial 
t-derivative by a "dot"): 

(2.1) P + P' z z = 0. 

Since ip(t) maps T t onto the unit circle, we have z(t) = $?(t,w), where 
\&(t, w) = l ^>~ 1 (t, z). We also have on T t , by differentiating P(z(w),t) = 
w n with respect to w, 

(2.2) P' z ■ z w = nw 11 - 1 
or 

(2.3) wz v 



or 

nw n nP 
P' ~ ~P>~' 



z 



From ( 12. ip , conjugating, we infer 

P 

(2.4) z = - = . 
v ' pi 

z 

Parametrize the unit circle by w = e t0 , < 9 < 2tt. Then, from (12. 3p . 
it follows that we have on T t (since (z(t) = z(w,t) = z(w(9,t)))), 

. 1 dz nP 

(2.5) r e:= m = ZwW = T>- 

Combining (12 .4p and (12.51) yields (3? stands for the real part): 

(2.6) 3? z e 
Also, 

(2.7) B (l^j- ) = 9 (Izo) = 1, 
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where in the last equality we used the hypothesis that the lemniscates 
r t := {z(t, 9)} satisfy the main equation (ll.2p of Laplacian growth 
processes— cf. [8, §4]. Hence, (J22D, (|Z3D and imply that 

-1 

PL ~' u I " \ PIP 7 I ~ '~n' 




z 




(2.8) 

Or, we can rewrite (12. 8p as 

(2.9) ft(PF) = -^\ p 'z\ 2 - 

Thus, we are finally arriving at 

(»o) |(|ph = _±|^. 

Therefore, (I2.1UI) holds on the lemniscates T t = {\P(z,t)\ = 1} that 
are assumed to be interfaces of a Laplacian growth process. Now the 
theorem follows from the following. 

n 

Lemma 2.1. Let t be the time variable, P(z,t) = a(t) Y\{z — Xi{t)) , 

i 

be a "flow" of n-degree polynomials. Assume that the lemniscates 
T t := {\P{z, t)\ = 1} all have connected interiors {\P(z,t)\ < 1} and a 
generalized equation (I2.10p holds on T t ; i.e., 

(2.11) j t {\P{z^)-c{t)\P' z {z,t)\ 2 = ^ 

where the function c(t) is real-valued, depends on t only and, hence, is 
a constant on T t . Then, n = 1, Ai = \i(t) = const and T t is a family 
of concentric circles centered at X±. 

Proof of the Lemma. Our hypothesis implies that all polynomials 
\P(z,t)\ 2 — 1 are irreducible. Hence, using Hilbert's Nullstellensatz 
(e.g., cf. [2], Proposition 3.3.2), we infer from (12. lip that 

(2.12) j t (|P(z,t)| 2 ) - c(t) \P'{z,t)\ 2 = B(t) (|P(M)| 2 - 1) • 

Equation (12.121) holds for all z G C and for an interval of time t, and for 
each t, both sides are real-analytic functions in z and z. Hence, we can 
"polarize" ( I2.12p . i.e., replace z by an independent complex variable 
£. (This is due to a simple observation: real-analytic functions of two 
variables are nothing else but restrictions of holomorphic functions in 
z, ^-variables to the plane {£ = z}. Hence, if two real-analytic functions 
coincide on that plane, they coincide in C 2 as well.) Denoting by 
P* the polynomial whose coefficients are obtained from P by complex 
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conjugation, we have (12.121) in a "polarized" form holding for G 

(P2. 



(2.13) - (P(z,t)P*(U)) -c(t) (P>(z,t) ■ (P%(U)] 

= B(t) (P(z,t)P*(Z,t)-l). 

Now let us denote by kj the multiplicity of the root Xj (t) of the polyno- 
mial P(z,t), so that there are m < n distinct roots and Y^jLi % = n - 
Since 



P(z } t) = a(t)l[(z - X 3 (t)) k 
1 

m 

P*(U) = a(t)H(t-W)) 
1 

dividing by P(z,t)P#(£,i) we obtain: 



(2.14) 2SR[-]-V fcA(t) fc^L 



-c(t) 



V — 



B(t) (l 



P(z,t)P#(U)J' 



Integrating ( 12.141) along a small circle centered at Xj(t), so that it does 
not enclose other zeros of P, yields for all £: 



(2.10) kjXj(t)- , 



where 



i 



(fc-r-i)! 



P(z,t) 



. Letting £ -> oo in (12451) 



implies that Aj(t) = for all j = 1, . . . , n. In other words, the "nodes" 
Xj(t) of all the lemniscates T t are fixed, i.e. do not move with time. 
So, 



(2.16) 



P(z, t) = a(t) Y[(z- Xj) = a(t) Q{z) 
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Substituting ( 12. 16ft into ( 12. 13ft . we obtain 

(2.17) j t (|a| 2 ) Q(z)Q#(0 - c(t)\a\ 2 Q' z (Q% 

= B(t)(\a\ 2 Q(z)Q#(0-l). 

Comparing the leading terms (i.e., the coefficients at z n £ n ) in ( 12.171) 
yields 

(2.18) j t (|a| 2 ) = B(t) \a\ 2 . 
Therefore, 

(2.19) c(t) \a\ 2 Q' z (Q#)[ = B(t), 

and thus degQ' z = 0, i.e., n = degP = 1. The proofs of the Lemma 
and the Theorem are now complete. □ 

3. Extending the theorem to growth processes invariant 

under time reversal 

First, let us note that any boundary T t is an equipotential line of the 
logarithmic potential 

(3.1) Hz) =\og\P n (z,X i (t))\ 2 - 

The boundary velocity of the general growth process defined in ( 11.31) 
can now be expressed as V(z) = x(z)V$, z G r f , x{ z ) £ 

As indicated in the Introduction, invariance under time-reversal is 
defined here in the following way: if the boundary T t+ dt is the image of 
r t under a map f {t>dt) : z t G V t i-> z t+dt G T t+dt , then f( t +dt-dt)°f(t,dt) = 
I. That means that the normal at z t+ dt G T t+dt must be parallel to the 
normal at z t G T t , which shows that T t +dt is perpendicular at every 
point to gradient lines of and is therefore a level line of $. The 
displacement of the point z t becomes 

zt+dt ~ Z-t = x(z)V^(z t )dt. 

Denoting by 



E = V$ = 2<9$ = 2 • P "( z > Xi( ty[ 

Pn(z,Xi(t)) 

the gradient of the logarithmic potential and by r = zt, conservation 
of the normal (or gradient) direction becomes 

E{f+ xE{r)dt) = fi{z)E(f), 
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where //(z) = 1 + m(z)dt, m(z) = 0(1), m(z) 
in the infinitesimal time interval dt, 



so after expanding 



(E ■ V)£(f) 



771(2) 



E(r) 



Remark 3.1. The proportionality relation indicated above carries also 
the following physical significance: the dynamical system that we study 
is of frictional type, where the acceleration field (proportional to the 
force, or gradient of Green's function) is also proportional to the veloc- 
ity. In other words, the transport derivative (or Lie derivative) of the 
velocity field must be parallel to the velocity itself: 

^y=[ V od-do i 9 ]V = (V ■ V)V = x\x(E ■ V)E + (E ■ Vx)E] 

is parallel to V and therefore, to E. 

Ed + Ed, we 



In complex notation, using the fact that (E ■ V) 
obtain 



S(z 



P n (z,\i(t)) v J P n (z,\i(t)) \P n (z,\i(t)) 
which (after multiplying both sides by E(z)) reduces to 



5(z)e 



or 



(3.2) 



P n (z,Xi(t))_ 

Pn(z, A*(t)) 



0, W)z E T t . 



[P^Ht))\ 

We note that, since E(z) = 1P' n jP n is the gradient of the Green's 
function for Q t and fl t is simply connected, it cannot vanish anywhere 
in fl t U F t , so all the zeros of Pn(z), denoted by k = 1, 2, . . . , n — 1, 



are found inside the domain Vt' t . Then 



Pnjz, \j(t)) 



1 n— 1 

1 + E 

n ^ 



k=l 



\2- 



with Ak constants. The imaginary part of this expression coincides 
with the imaginary part of an analytic function in Q t , that is bounded 
there, so the condition (13. 2p can only be satisfied if the function is a 
constant. Since at z — > oo it vanishes, it follows that 



Pn{z, \i(t)) 



1 

J 

n 



which means that boundaries T t can only be concentric circles. 
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4. Concluding Remarks 

(1) It is plausible that the result can be extended to rational lemnis- 
cates F t := {\R (z, t)\ = 1}, where R (z, t) are rational functions 
of degree n where all the zeros are inside T, while all poles are 
in the unbounded component of C \ T t . 

(2) It is well-known that arbitrary "shapes" , i.e. Jordan curves can 
be arbitrarily close approximated by both lemniscates (Hilbert's 
theorem - cf. [14]) and quadrature domains [1]. At the same 
time our results imply that there are fundamental differences 
between these two classes of curves. We think it is interesting 
to pursue these observations in greater depth. 

(3) From the argument in §|3]we can extract more. Suppose a family 
of Jordan curves, {r t } t>0 evolves by the flow along the velocity 
field V(z) according to (11.31) . Assuming the invariance under 
time-reversal, the argument of §|3]can be used to prove that x — 
const, i.e. the process is that of Laplacian growth. Invoking 
now well-known results on standard Hele-Shaw flows, we can 
at once conclude, e.g., that the process (11.31) continues for all 
times t > 0, i.e., the curves {I\} move out to infinity such that 
Ut>t o r t = C \ Q to , if an d only if the initial curve r is an ellipse 
and all the curves {T t } are also ellipses homotetic with r - 
cf. [3], also cf. [6]. 
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